A probabilistic model for the numerical solution of initial value problems

نویسندگان

  • Michael Schober
  • Simo Särkkä
  • Philipp Hennig
چکیده

Like many numerical methods, solvers for initial value problems (IVPs) on ordinary differential equations estimate an analytically intractable quantity, using the results of tractable computations as inputs. This structure is closely connected to the notion of inference on latent variables in statistics. We describe a class of algorithms that formulate the solution to an IVP as inference on a latent path that is a draw from a Gaussian process probability measure (or equivalently, the solution of a linear stochastic differential equation). We then show that certain members of this class are connected precisely to generalized linear methods for ODEs, a number of Runge–Kutta methods, and Nordsieck methods. This probabilistic formulation of classic methods is valuable in two ways: analytically, it highlights implicit prior assumptions favoring certain approximate solutions to the IVP over others, and gives a precise meaning to the old observation that these methods act like filters. Practically, it endows the classic solvers with ‘docking points’ for notions of uncertainty and prior information about the initial value, the value of the ODE itself, and the solution of the problem.

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Implicit One-step L-stable Generalized Hybrid Methods for the Numerical Solution of First Order Initial Value problems

In this paper, we introduce the new class of implicit L-stable generalized hybrid methods for the numerical solution of first order initial value problems. We generalize the hybrid methods with utilize ynv directly in the right hand side of classical hybrid methods. The numerical experimentation showed that our method is considerably more efficient compared to well known methods used for the n...

متن کامل

The symmetric two-step P-stable nonlinear predictor-corrector‎ ‎methods for the numerical solution of second order‎ ‎initial value problems

In this paper‎, ‎we propose a modification of the second order method‎ ‎introduced in [‎‎Q. Li and ‎X‎. ‎Y. ‎Wu‎, A two-step explicit $P$-stable method for solving second order initial value problems‎, ‎textit{‎Appl‎. ‎Math‎. ‎Comput‎.}‎ {‎138}‎ (2003)‎, no. 2-3, ‎435--442‎] for the numerical solution of‎ ‎IVPs for second order ODEs‎. ‎The numerical results obtained by the‎ ‎new method for some...

متن کامل

Trigonometrically fitted two-step obrechkoff methods for the numerical solution of periodic initial value problems

In this paper, we present a new two-step trigonometrically fitted symmetric Obrechkoff method. The method is based on the symmetric two-step Obrechkoff method, with eighth algebraic order, high phase-lag order and is constructed to solve IVPs with periodic solutions such as orbital problems. We compare the new method to some recently constructed optimized methods from the literature. The numeri...

متن کامل

NON-STANDARD FINITE DIFFERENCE METHOD FOR NUMERICAL SOLUTION OF SECOND ORDER LINEAR FREDHOLM INTEGRO-DIFFERENTIAL EQUATIONS

In this article we have considered a non-standard finite difference method for the solution of second order  Fredholm integro differential equation type initial value problems. The non-standard finite difference method and the composite trapezoidal quadrature method is used to transform the Fredholm integro-differential equation into a system of equations. We have also developed a numerical met...

متن کامل

‎A Consistent and Accurate Numerical Method for Approximate Numerical Solution of Two Point Boundary Value Problems

In this article we have proposed an accurate finite difference method for approximate numerical solution of second order boundary value problem with Dirichlet boundary conditions. There are numerous numerical methods for solving these boundary value problems. Some these methods are more efficient and accurate than others with some advantages and disadvantages. The results in experiment on model...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

عنوان ژورنال:
  • CoRR

دوره abs/1610.05261  شماره 

صفحات  -

تاریخ انتشار 2016